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Abstract 

The quantification of ttie stiffness of tubular biological structures is often obtained, both 
in vivo and in vitro, as the slope of total transmural hoop stress plotted against hoop strain. 
Total hoop stress is typically estimated using the "Laplace law." We show that this procedure 
is fundamentally flawed for two reasons: Firstly, the Laplace law predicts total stress incor- 
rectly for biological vessels. Furthermore, because muscle and other biological tissue are closely 
volume-preserving, quantifications of clastic modulus require the removal of the contribution to 
total stress from incompressibility. We show that this hydrostatic contribution to total stress 
has a strong material-dependent nonlinear response to deformation that is difficult to predict 
or measure. To address this difficulty, we propose a new practical method to estimate a me- 
chanically viable modulus of elasticity that can be applied both in vivo and in vitro using the 
same measurements as current methods, with care taken to record the reference state. To be 
insensitive to incompressibility, our method is based on shear stress rather than hoop stress, 
and provides a true measure of the elastic response without application of the Laplace law. We 
demonstrate the accuracy of our method using a mathematical model of tube inflation with 
multiple constitutive models. We also re- analyze an in vivo study from the gastro-intestinal lit- 
erature that applied the standard approach and concluded that a drug-induced change in elastic 
modulus depended on the protocol used to distend the esophageal lumen. Our new method 
removes this protocol-dependent inconsistency in the previous result. 

1 Introduction 

Although originally developed to estimate mathematically the surface tension at the air-liquid 
interface in bubbles as a function of geometry and bubble pressure, the "Laplace law" is used widely 

'Submitted for publication in Mathematical Medicine & Biology (a journal of the IMA) on 4 February 2013, 
manuscript ID: MMB-13-005. 
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in the medical community to estimate the average hoop stress in the wah of biological vessels such as 
blood vessels and the muscle-driven conduits of the gastrointestinal (GI) tract (Basford, 2002). In 
the GI tract, for example, the Laplace law is commonly used to estimate average hoop stress across 
the esophageal wall* as a function of distension and manometrically measured pressure within the 
esophageal lumen (Biancani et al., 1975; Fr0kja3r et al., 2006; Ghosh et al., 2005; Gregersen et al., 
2002; Pedersen et al., 2005; Takeda et al., 2003, 2002, 2004; Yang et al., 2004, 2006a,b,c, 2007). 
The Laplace law is also commonly applied in the cardiovascular system (see, e.g., Kassab, 2006). 

In addition, to estimate material stiffness, the wall-average hoop stress of the vessel wall is 
commonly plotted against a stretch measure and the slope of the cTirve interpreted as a modulus 
of elasticity or wall stiffness.^ For example, Takeda et al. (2003, 2004) determine a modulus of 
elasticity by plotting the circumferential second Piola-Kirchhoff stress against the circumferential 
component of the Green strain,^ and Gregersen et al. (2002) plot the Cauchy circumferential stress 
against a circumferential stretch ratio, where in both cases stress is estimated using the Laplace law. 
A major advantage of this approach is that, with the Laplace law, hoop stress and deformation can 
be estimated, in vivo as well as in vitro, from measurements of intraluminal pressure and luminal 
cross sectional area with concurrent manometry and endoluminal ultrasound (Gregersen et al., 
2002; Takeda et al., 2003, 2002, 2004; Yang et al., 2004). With this approach, luminal wall stiffness 
is estimated and compared across different clinical groups such as patient /healthy, old/young and 
pre/post therapy, or over time. 

We show that this approach to measuring a modulus of elasticity is flawed for two reasons. 
Firstly, we show that the Laplace law is invalid and in serious error as an estimate of total hoop 
stress when applied to biological vessels with physiological deformations. Secondly, we demonstrate 
that the slope of average total hoop stress plotted against an appropriate strain measure does not 
represent a modulus of elasticity, both in principle and in practice, because, to measure stiffness, 
the hydrostatic component of the stress due to incompressibility must be removed. Clearly, there 
remains a great need for a practical approach to estimate average material stiffness of biological 
vessels in vivo and in vitro. We propose an alternative strategy with which a true modulus of 
elasticity can be estimated from data obtained via current in vivo and in vitro experimental methods 
and we demonstrate the validity of our new method using parameters for the esophagus. 

2 The Laplace law as an estimate of average hoop stress in bio- 
logical vessels 

Wall stiffness of biological vessels is commonly estimated by interpreting the slope of a total hoop 
stress against strain as an elastic modulus. Both in vivo and in vitro, the raw data used to generate 
the curve consist of measurements of the transmural pressure difference and of the inner and outer 
radii of the vessel. The Laplace law is then used to translate these measurements into the needed 
wall-average total hoop stress. In this study we address the following questions: 

1. Since the Laplace law is an approximation valid in the limit of thin- walled vessels, how well 
does the Laplace law approximate the total average hoop stress in biological vessels? 

2. Since muscle and other biological tissue are nearly incompressible and the total hoop stress 

*The esophageal wall is composed of circular and longitudinal muscle layers, and inner mucosal layers of mostly 
loosely connected tissue. 

^Hcrc the word "stiffness" is used, as in the medical literature, to mean "material resistance to deformation." 
*This strain is also called the Green-Saint Venant strain(Gurtin et al., 2010), the Green-Lagrange strain (Holzapfel, 
2000), or the Lagrangian strain (Bowen, 1989; Taber, 2004). 
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has both an elastic component and a component due to volume conservation (due to incom- 
pressibility), does the total average hoop stress (plotted against strain) approximate well the 
elastic response of the vessel (relative to the purely elastic component)? 

We address these questions by first reviewing the derivation of the Laplace law; then comparing 
the average hoop stress estimates obtained via the Laplace law with exact estimates based on two 
rigorous solutions. 



2.1 The Laplace law 



The Laplace law represents an equilibrium force balance across the wall of a tubular structure with 
an infinitesimally thin wall. Consider in Fig. 1 a circular cylinder that is distended in response to 
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Figure 1: A circular cylinder and the coordinate systems used to describe the reference and deformed 
configurations. 

a transmural pressure difference = Pi — Po across the wall, where pi and po are the inner and 
outer pressures, respectively. The cylinder's inner and outer radii are and Tq, respectively, so 
that the wall's thickness is i = — rj. If the geometry of the cylinder is axisymmetric and the 
material is isotropic, the Cauchy hoop stress ctoq depends only on the radial coordinate r and not 
on the axial coordinate z or azimuthal angle 9. Denoting by (•) the average of a quantity across 
the cylinder wall, the average hoop stress is 



{(yee) 



1 



aee dr. 



(1) 



Cutting the cylinder longitudinally as illustrated in Fig. 2, and balancing the resultant force per 
unit length due to pi over half the cylinder in the direction perpendicular to the cut (with magnitude 
2piri) with the corresponding resultant force due to po and {(Tqq) (with magnitude 2poro + 2{agg)t) 
one obtains: 

Ap 

Wi 

The Laplace law is the limit of Eq. (2) for t/ri — )• 0. For bounded values of po and pi, the term 
Ap/{t/ri) becomes dominant relative to po and the equilibrium force balance for a "thin walled" 
cylinder, the "Laplace law," becomes 
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Figure 2: Free-body diagram of half of a cross section of the pressurized cyhnder. The average 
force distribution acting on the cut is shown. 



The thin wall approximation, typically expressed as t/vi <C 1, could be naively interpreted in 
practice by taking t/vi < 0.1. However, the passage from Eq. (2) to Eq. (3) requires that po <C 
Ap/{t/ri), leading to the following requirement for a proper use of the Laplace law: 

* '^P ox 
- < — . 4 

U Po 

This condition is rarely met in biological vessels. For example, radial distention of the esophagus 
during the transport of food boluses between the pharynx and stomach are driven by pressure 
differences up to a maximum of ~ 25-30mmHg (Ghosh et al., 2005). Therefore, since po is close to 
atmospheric pressure, for the Laplace law to be applicable in the esophagus, the relative thickness of 
the esophageal wall must be extremely small: t/rj <C 0.04. This condition is almost never met in the 
human body. For example, extensive measurements of esophageal cross sections using endoluminal 
ultrasound before, during, and after esophageal distension with bolus transport (Nicosia et al., 
2001; Schiffner, 2004; Ulerich et al., 2003) indicate that t/ri typically ranges between a maximum 
value (t/ri)max ~ 0.7, in the pre-swallow resting state, to a minimum value {t/ri)iam ~ 0.1 in the 
most distended state. Thus the condition t/ri ^ 0.04 is never met, not even approximately. 

Analysis in blood vessels leads to a similar conclusion. The transmural pressure difference Ap, 
for example, is approximately 100 mmHg or less, depending on the vessel and degree of hypertension 
(Heagerty et al., 1993; Kim et al., 2002; Pries et al., 2005), so that Ap/po ^ 0.13 or less. Thus, 
t/ri must be much smaller than 0.13 for the Laplace law to be a reasonable approximation. Kim 
et al. (2002), using MRI, report normal t/r^ ratios in the coronary arteries around 0.5 in humans 
and much higher with hypertension. Pries et al. (2005) report values of t/ri from the literature for 
arterioles in rats ~ 0.21-0.26, while Heagerty et al. (1993) report in vitro values of normal small 
arteries in humans ~ 0.10-0.13, and for hypertensive small arteries ~ 0.16-0.20. None of these 
values come close to the requirement that t/ri be <^ 0.13 for the Laplace law to be an accurate 
approximation of total hoop stress. 

With the above in mind, we wish to provide a quantitative assessment of the accuracy of the 
Laplace law in incompressible materials. From this assessment, we can then present an objective 
critique of the practice of measuring a modulus of elasticity from a plot of such average stress 
against strain. This goal can be readily achieved by determining an exact solution of the inflation 
and extension problem for an elastic tubular vessel. Here, we considered two that are particularly 
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simple: the first for an isotropic nco-Hookean material, and the second for a simplified Mooney- 
Rivlin material. In both cases, the constitutive model is chosen so to have a single constitutive 
parameter that could be estimated using the information measured in vivo. Whether or not the 
chosen models are appropriate for the modeling of biological tissue is a secondary concern here. 
The analysis presented here is meant to show in a quantitative way that current practices for the 
determination of tissue stiffness are inadequate and that the development of more rigorous methods 
for the task is justified. 

2.2 Inflation and Extension of Tubuleir Vessels: Exact Solutions for Two In- 
compressible Isotropic Materials 

To evaluate the accuracy of the Laplace law for biological vessels, we derived exact expressions 
for the stress field in an incompressible cylindrical vessel subject to prescribed axial stretch 
and transmural pressure difference Ap. Two different models were chosen to show that the same 
conclusions are reached in both cases. These models were also chosen so as to be able to derive 
solutions in closed form and no claims are made about their ability to represent the behavior of a 
particular biological tissue. 

Solutions to the infiation and extension of a circular cylinder has been discussed by various 
authors (see, for example, Humphrey, 2010; Ogden, 1997; see also the analysis by De Pascalis, 
2010). However, the present authors are not aware of solutions that show in detail the full stress 
field in the form needed for the analysis at hand. For this reason and for the sake of completeness, 
the relevant details of our solution are reported in the Appendix. 

Kinematics. Consider the inflation and extension of a circular cylinder with uniform wall thick- 
ness (Fig. 1). Points in the reference configuration are identified by the cylindrical coordinate 
coordinates {R, @, Z) (the origin is on the longitudinal axis of the cylinder) ranging as follows: 

Ri<R<Ro, < G < 27r, and < Z < L, (5) 

where Ri, Rq, and L are the internal radius, external radius, and length of the cylinder, respectively. 
As discussed later in detail, we assume that the reference configuration is unloaded in the sense 
that the pressure difference across the cylinder's wall is equal to zero. Points in the deformed 
configuration are identified by the coordinates {r,9,z). By inflation and extension we mean a 
deformation such that 

r = r{R), 9 = e, z = \zZ, (6) 

where r{R) is at least twice differentiable. The quantity > 0, uniform across the wall, is the 
prescribed longitudinal stretch of the cylinder. We denote by ri and rg the values of r for i? = Ri 
and R = respectively. We use two sets of orthonormal base vectors, namely (m_r, uq, uz) and 
u,f,UQ,Uz) , tangent to the corresponding coordinate lines of the {R,Q,Z) and {r,6,z) coordinate 
systems, respectively. When the material is incompressible the volume is preserved so that we must 
have 

r'^A, = l ^ X4r^-r^)=R^-Rl (7) 

where = d • /dR, and where r = for i? = Prom the second of Eqs. (7) we see that, when 
Xz = 1 (i.e., no longitudinal shortening or extension), the cross-sectional area of of the cylinder is 
also preserved during deformation. 
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For future reference, the deformation gradient F, F ^, the left Cauchy-Green strain tensors B, 
and its inverse B~^ axe: 

F = Ur'SlUR+ —U0'S)Ue + XzUz'Sl Uz, (8) 

Xz r R 

r R 1 

F~-^ = A2— (g) H — lie ^ iie + ® iiz, (9) 
R r Xz 

1 i?2 j,2 

B = Ur<^ Ur + j^ue '^Uo + XlUz<^ Uz, (10) 

B ^ = Xl^Ur<^Ur+ ^Ue'^U0 + ^Uz<^Uz- (H) 

Stress field. To determine the stress field in the deformed state it is necessary to define the 
material through constitutive equations. As mentioned earlier, our choice of constitutive equations 
is motivated by our desire to quantify the accuracy of the Laplace law as an estimate for the total 
average hope stress. We choose two constitutive relations from which full and exact stress solutions 
can be obtained: a neo-Hookean model and a simplified Mooncy-Rivlin model. By modeling the 
esophagus, for example, as a homogeneous, isotropic, and incompressible hyperelastic material, 
we neglect rate dependence (Yang et al., 2006c), the layered structure of the esophagus (Ghosh 
et al., 2008; Nicosia and Brasseur, 2002; Yang et al., 2006a,b, 2007), and the potential existence of 
heterogeneity and non-hydrostatic residual stress fields. 

Since average intrathoracic pressure is approximately atmospheric, and because the esophageal 
reference state has been taken as the unloaded, the inner and outer pressures pi and Po are approx- 
imately equal, in the reference state we set pi and po equal to atmospheric pressure. In addition, 
in our analysis we consider equilibrium states under negligible body forces. Because the material is 
incompressible, if all of the material constitutive parameters are known and the (volume-preserving) 
deformation prescribed, the stress solution is completely known to within a Lagrange multiplier q, 
which is found by enforcing static equilibrium and boundary conditions. In addition, limiting our 
attention to incompressible isotropic materials with undistorted reference configurations* (Bowen, 
1989), the only nontrivial equilibrium equation of the problem, written in terms of the Cauchy 
stress, has the following form: 



dr 



+ "^"^ = 0. (12) 



Neo-Hookean materiaL The Cauchy stress constitutive response function for a homogeneous 
isotropic neo-Hookean material with undistorted reference configuration is 

(T = -q\ + fi{B-\), (13) 

where q is the Lagrange multiplier required for the satisfaction of incompressibility. In our problem, 
B is given by Eq. (10). Since n is constant, we absorb it into q and rewrite Eq. (13) with the 
substitution q = q + ji: 

cr = -g|+^B. (14) 

Enforcing static equilibrium and boundary conditions, we find the following expression for q (see 
Eq. (A. 6) in the Appendix): 



(g)NH =Po + 1^ 



1 fRl , , 1 , frRo 

H FT + ^ In 



2A? \ / A^ \ r„R 



(15) 



*The main implication of this assumption is that if a residual stress field exists, this field must be hydrostatic. 
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where the subscript 'NH' stands for 'nco-Hookcan.' The quantity ^ can be expressed in terms of 
quantities that are measured and/or controlled in an experiment (see Eq. (A. 8) in the Appendix): 



Ap 



with 



j{ri,ro,Ri,Ro,X 
1 
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2A2 V 
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^ 1 iJ^jRo 

K \roRi 



(16) 



(17) 



Combining Eqs. (10) and (15) into Eq. (14), we have that the hoop stress a^Q for the neo-Hookean 
material is given by (see Eq. (A. 9) in the Appendix) 



NH 



-Po - 



1 



R„ R 



2A? \ r"^ ^ 



^ 1 , frRo 



(18) 



In practice, the empirical measurements of wall stiffness are made using wall-average quantities. 
This is precisely the case when the Laplace law is used since this law provides an estimate of 
average of the hoop stress across a vessel's wall. For this reason, we now use the pointwise results 
in Eqs. (15) and (18) to obtain their average across the wall of the cylinder. Specifically, letting 

p = ^rf — R^/\z = i/r^ — R^/Xz, and using Eqs. (A.20)-(A.25) in the Appendix and simplifying, 
we have 



, , p R^ iJ, 



+ 



(<76»e)NH 



Po + 



■ n) 



pin 
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TnRi 



R? . , fviRo 

ToRi 



(19) 
(20) 



Again, it is important to keep in mind that, when performing an experiment in which the quantities 
Ap, Ro, Tj, To and Xz are measured and/or controlled, the quantity p in Eqs. (19) and (20) can 
be determined using Eqs. (16) and (17). 

Mooney-Rivlin material. The Cauchy stress for an incompressible homogeneous isotropic 
Mooney-Rivlin material with undistortcd reference configuration is 



-gl + /xi(B-l)+/i2(l-B-i), 



(21) 



where q is a Lagrange multiplier for the enforcement of incompressibility, and where p\ and p2 are 
constants. With the substitution q = p\ — P2, Eq. (21) becomes 



(T = —q\ + piB — /X2B ^ 



(22) 



In our problem, the B and B~^ are given by Eqs. (10) and (11), respectively. For pi = p and p2 = 0, 
Eq. (22) yields the neo-Hookean model considered above. To reduce Eq. (22) to a one-parameter 
model different from the neo-Hookean, we therefore consider a "reduced Mooney-Rivlin" model 
with pi = 0: 



cr = -q\ + 2p2£, 



(23) 
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where 



g + /i2 and e = 1(1 — B 



(24) 



The tensor e is often referred to as the Eulerian strain tensor (c.f. Taber, 2004) although other 
names used are the Hamel (c.f. Tabcr, 2004), Almansi (c.f. Taber, 2004), and Euler-Almansi (c.f. 

Holzapfcl, 2000) strain tensors. Since B^^ is known (see Eq. (11)), the stress cr in Eq. (22) is 
completely known once the field q is known. The latter can be shown to be equal to (see Eq. (A. 15) 
in the Appendix) 



(9)rmr = Po + /^2A^ 



1 

2A! 



+ 



1 , 

+ T-ln 



rRo 
roR 



W ^ a! 72" 



(25) 



where the subscript 'RMR' stands for 'reduced Mooney-Rivlin'. Similarly to the neo-Hookean case, 
II2 can be expressed in terms of quantities that are measured and/or controlled in an experiment. 
Specifically, we have (see Eq. (A. 17) in the Appendix) 



Ap 



>^h{ri,ro,Ri,RoAz)' 



(26) 



where 7 is given by Eq. (17). Substituting m = and Eq. (25) into Eq. (22), the hoop components 
of the Cauchy stress tensor is 



RMR 



-Po - M2A 
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2A! 



it„ it 

— H 

J.2 J.2 



+ t'-Iu 



rRo 
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(27) 



Proceeding as shown for the neo-Hoookean material, the transmural averages of (g)RMR) and 
(o-6>6i)rmr are 



/ \ I /^2 , /^sA^ 

(g)RMR = Po + y ^ + 



P 
Ton 



1^2 Rq , fJ'2^z 



r„ri 



(fee) RMR 

where we recall that p = J rf — R^ / \z = -^r^ — R^/Xz- 



+ ^ln^ + l 



ToRi 



(28) 
(29) 



3 Inadequacy of the Laplace law under physiological conditions 

Using the stress solutions for the two elementary models considered, we now quantitatively assess 
the accuracy of the Laplace law under normal physiological conditions in the human esophagus. 
Typical cross-sectional areas and thicknesses of the muscle layers (muscularis propria) of the esoph- 
agus have been measured using intraluminal ultrasound (Miller ct al., 1995; Mittal et al., 2006; 
Nicosia and Brasseur, 2002; Schiffner, 2004). It is well known that the intrathoracic pressure os- 
cillates with respiration by several mmHg and with a mean value of a couple of mmHg below 
atmospheric pressure (Christensen, 1987). A value for the constant 112 was estimated by Schiffner 
(2004) to be approximately 60 mmHg. This value was used by Nicosia and Brasseur (2002) and 
Ghosh et al. (2008) to model esophageal muscle. In the current analysis we use the following pa- 
rameters as representative of the esophagus in the reference state (when pi = po), before infiation 
by a swallowed bolus of fluid: 

Ri = 2.6mm, Ro = 4.4mm, /x = /X2 = 60 mmHg, po = 760 mmHg, T/Ri = 0.692, (30) 
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where T/Ri is the muscle layer thickness ratio in the reference state. Therefore, we evaluate the 
Laplace law in the range 

< t/n < 0.7. (31) 

As alluded to in the Introduction, for the Laplace law to be accurate to within 10%, it should 
only be applied when the thickness ratio t/ri is approximately between and O.lAp/po- Although 
one normally swallows liquids and solids as boluses less than 10 ml in volume, clinical studies are 
commonly carried out with bolus volumes up to 20 ml. Fluoroscopic and endoluminal ultrasound 
measurements of luminal distension of the esophageal body with bolus volumes up to 20 ml (Chris- 
tensen, 1987; Ghosh et al., 2005; Nicosia and Brasseur, 2002; Shi et al., 2002) indicate a lower 
bound for t/ri of about 0.12. Furthermore, it is well known that longitudinal muscle fibers in the 
esophageal wall contract during bolus transport (Dodds et al., 1973) creating both local and global 
longitudinal shortening of the esophagus in the range 0.33 < < 1.1 (Dai et al., 2006; Mittal 
et al., 2006; Nicosia and Brasseur, 2002; Shi et al., 2002). Therefore, the range of thickness ratio, 

0.12 < t/ri < 0.69, (32) 

will be referred to as the physiological range. While we do not focus on longitudinal shortening 
in this study, to observe how may affect the response of the system we present results for 
A, = 0.8,1,1.1. 



Stress Estimates. In Fig. 3 we plot Ap/po as a function of the nondimensional thickness t/ri for 
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Figure 3: Ap/po vs. wall thickness ratio t/ri. The behavior of the NH and RMR models is shown 
via solid and dashed lines, respectively. Three values of A^ are plotted: 0.8 (red), 1 (blue), and 
1.1 (green). The NH and RMR lines for = 1 coincide. The physiological range is indicated by 
vertical blue lines. The vertical grey strip on the left represents the "Laplace law" range defined as 
< t/ri < O.lAp/po ^ 0.004. 

the NH and RMR models generated using Eqs. (16), (17), and (26). Subject to the constraint in 
Eq. (4) , we expect the Laplace law to provide a reasonable estimate of the average hoop stress when 
< t/ri < O.lAp/po. Fi gure 3 shows that, for the NH and RMR models, Ap/pQ varies between 
zero and approximately 0.04, which implies that the Laplace law estimates of the average hoop 
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Figure 4: Ratio of exact transmural average stress {(Jqq) to the corresponding average (o-6»e)LL 
approximated by the Laplace law (Eq. (33)). The behavior of the NH and RMR models is shown 
via solid and dashed lines, respectively. Three values of are plotted: 0.8 (red), 1 (blue), and 1.1 
(green). The NH and RMR lines for = 1 coincide. The physiological range is to the right of the 
vertical cyan line. The vertical grey strip on the left represents the "Laplace law" range defined as 
< t/n < O.lAp/po 0.004. 

to the Laplace law approximation {age) ll in Eq. (3) is plotted against the thickness ratio t/ri. For 
the NH and RMR models, we have 

(am) (l-^ for the NH model, 
{(^ee)LL I 1 for the RMR model. 

Figure 4 shows that the the Laplace law provides accurate estimates of the average hoop stress 
only when t/ri is less than about 0.001. Not only are the stress estimates in the physiological range 
badly in error, they are of the wrong sign. This is as a result of neglecting po in the full force 
balance, Eq. (2). That is, the full force balance and the Laplace law differ by the constant po, 
which is large relative to {(Tgg)i,i,, Eq. (3). 

Since {crge)LL and (aog) differ by the constant po, and since po must be estimated in any event to 
evaluate Ap = Pi — Po from a measurement of pi , one could argue that the error in the application 
of the Laplace law to biological vessels is easily rectified by using the full force balance, Eq. (2). 
However we show next that the characterization of the elastic response of an incompressible material 
such as muscle using total stress is incorrect, both in principle and in practice. 

4 Pollution of true elastic response by incompressibility 

In modeling the elastic response of incompressible materials, the total stress should be decomposed 
as 

o■ = -gl + o■^ (34) 
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where q is the Lagrange multipHer for the enforcement of incompressibihty and cr*^ is the part of 
the stress response that is governed by deformation (i.e., the part that is known when an appro- 
priate measure of deformation is known). To characterize the elastic response of the material, one 
must characterize the behavior of cr*^ rather than cr. Since {(Joe), the transmural average of the 
hoop component of cr, is much easier to determine experimentally than the (cqq), we quantify the 
difference between (ago) and {cTgg) (by quantifying (q)), and ask whether the Laplace law estimate 
(o'e6»)LL might approximate (cr^g), since we have established that it cannot be used to model {o-go). 
We use the solutions we have obtained for the NH and RMR models to address these issues. 

The expressions in Eqs. (19) and (28), nondimensionalized with respect to po, are plotted in 
Fig. 5. This figure shows that the transmural average of the Lagrange multiplier is strongly ma- 
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Figure 5: Transmural average of the Lagrange multiplier q vs. wall thickness ratio t/rj. The 
behavior of the NH and RMR models is shown via solid and dashed lines, respectively. Three 
values of are plotted: 0.8 (red), 1 (blue), and 1.1 (green). The physiological range is between 
the two vertical cyan lines. 

terial dependent, is comparable to the external pressure, and varies significantly with deformation 
especially when compared to the variation of the transmural pressure difference over the same range 
(cf. Fig. 3). Thus the effect of incompressibihty is a significant fraction of the transmural average 
of the total hoop stress and should not be ignored. This is shown explicitly in Fig. 6, where, over 
the entire physiological range, the quantity {q) is of order or greater than the average total hoop 
stress. In fact, Fig. 6 indicates that (agg) is mostly a reflection of {q)\ 

We conclude that (agg) is inappropriate for the determination of a modulus of elasticity. Con- 
sider, for example, the common empirical approach in which the transmural average of the total 
circumferential stress is plotted against an appropriate strain measure, and a modulus of elasticity 
is estimated from the slope of the curve, as shown in Fig. 7. Following the approach demonstrated 
by Takeda et al., 2003, 2004, in this figure, we plot {Sqq), the average of the circumferential com- 
ponent of the second Piola-Kirchhoff stress S (see Eqs. (A. 10), (A. 19), and Eqs. (A.20)-(A.25) in 
the Appendix) as a function of {Eqq), the corresponding average component of the Lagrangian 
strain E, where S and E are defined as follows:* 

S = (detF)F-VF-'^ and E=i(F'^F-l) (35) 
*The tensor E is the work conjugate of S. Also, using Eq. (8), we have that -Eee = [{r^ /R^) — l]/2- 
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Figure 6: Ratio of the transmural averages of q and agg vs. wall thickness ratio t/ri. The behavior 
of the NH and RMR models is shown via solid and dashed lines, respectively. Three values of 
are plotted: 0.8 (red), 1 (blue), and 1.1 (green). The physiological range is between the two vertical 
cyan lines. All plotted curves approach infinity as t/ri — )• 0. 



In addition, we show the contributions to {Sqq) due to the Lagrange multiplier q for the NH 
and RMR models. For the RMR model, the contribution due to the Lagrange multiplier is very 
hard to distinguish from (S'ee) so the slope in this case would primarily measure a "modulus of 
incompressibility," rather than a modulus of elasticity. In the NH case the situation is even worse. 
In fact, from Eq. (9), (10), the first of Eqs. (35), and (14), we have that (S'ee) = -{iR^/r^)q) + 
where — ((i?^/r^)g) is the contribution to (See) due to the Lagrange multiplier g, given by the 
solid red curve in Fig. 6. Since /i is constant, for the NH model, the slope of the (S'ee) curve is 
identical to that of the solid red line; i.e., it is purely a reflection of incompressibility. The point is 
that, for incompressible materials, the slope of (See) vs. its conjugate strain has very little to do 
with the elastic response of the material. The situation does not improve if other stress measures 
are used along with other corresponding strain measures. 

From this analysis, we conclude that the elastic behavior of an incompressible material cannot 
be characterized using the total stress due to the strong influence of the hydrostatic contribution. Is 
it possible, however, that the Laplace law estimate, which does not approximate the wall averaged 
total hoop stress, might provide a fortuitous approximation to the slope of the elastic component of 
the hoop stress {(Jgg) when plotted against an appropriate strain measure? For this to be the case, 
the ratio of (<T^g) to {crgg)i,i^ must be roughly constant and of order 1 during deformation by inflation. 
Figure 8 shows that this is not the case. Not only does the ratio of (c^g) to (o"eg)LL strongly vary 
during inflation of the vessel, its nonlinear variation is strongly material dependent. Thus, another 
approach is needed to replace the incorrect application of the Laplace law, an approach that "filters 
out" the contribution to total stress due to incompressibility (as well as residual stress) in order 
to quantify a true elastic modulus, or "stiffness" of tubular biological vessels. We present a new 
practical method to estimate a true modulus of elasticity using the same data obtained with the 
current approach based on the Laplace law. 
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Figure 7: Transmural average of the hoop component of the total second Piola-Kirchhoff stress for 
= 1 as a function of the wall averaged hoop component of the Lagrangian strain (blue). The 

results for the NH and RMR models are indicated by solid and dashed lines, respectively. For 
= 1 the NH and RMR (S'ee) curve coincide. The red lines show the contribution to (/See) due 

to the Lagrange multiplier. The physiological range is delimited by the two vertical cyan lines. 



5 A new approach to estimate a true modulus of elasticity for 
biological tubular vessels 

There is a great need for a method to estimate sensible stiffness measures of tubular muscle, 
in vivo as well as in vitro, that can be properly interpreted as moduli of elasticity in order to 
quantify changes associated with abnormality and disease. Although, in general, the structure 
of the luminal walls is not homogeneous or isotropic — in the GI tract, for example, muscle fibers 
are aligned azimuthally and longitudinally — for such clinical and therapeutic applications it is 
acceptable to use homogeneous isotropic models if they produce relatively simple measures that 
represent a true elastic response to deformation. A common current method applies an isotropic 
model and interprets muscle stiffness from the slope of wall-averaged total hoop stress, estimated 
from the Laplace law with pressure and geometry measurements, plotted against an appropriate 
circumferential strain measure. We have shown that this approach does not provide a sensible 
measure of an elasticity modulus as a quantification of muscle stiffness. 

There are two essential difficulties with the common approach. Firstly, the Laplace law does 
not approximate the true hoop stress for biological tissue and should be avoided in general. More 
importantly, to quantify the true elastic response of incompressible materials, the contribution to 
the stress due to incompressibility must be removed. If the influence of incompressibility were 
minor, in principle, the total stress plotted against an appropriate strain might retain a sufficient 
portion of the elastic response to serve as an approximate measure of stiffness. However, we have 
shown that incompressibility can be a dominant influence. 

In this section we propose a new approach to estimate muscle stiffness that is not influenced by 
incompressibility. The basic observation is that the shear response of an incompressible material is 
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Figure 8: Ratio of the transmural average of the hoop component of cr^ to the correspondent 
Laplace law approximation of the total stress. The behavior of the NH and RMR models is shown 
via solid and dashed lines, respectively. Three values of \z are plotted: 0.8 (red), 1 (blue), and 1.1 
(green). The physiological range is between the two vertical cyan lines. All of the curves pertaining 
to the NH and RMR models coincide; only the red curve is visible. 



unaffected by incompressibility. Thus, a true measure of elastic response can be obtained through 
an effective elastic shear modulus. The aim, therefore, is to develop a technique to estimate an 
effective elastic shear modulus within the standard experimental approach of inflation and extension 
of tubular organs in vivo as well as in vitro, i.e., using the same information that is used to estimate 
muscle stiffness as current methods. 

Our approach applies to incompressible tissue and assumes that the reference configuration 
is undistorted so that the residual stress is hydrostatic (Bowen, 1989) and can be absorbed into 
the Lagrange multiplier. To develop our approach, we relate simple shear and the inflation and 
extension test to define an effective shear modulus based on the same experimental method as 
current approaches. 

5.1 The simple shear test 

We review the simple shear test as a basis for our proposed approach to the determination of an 
effective shear modulus of tubular vessels from inflation and extension. The material presented 
here can be found in several textbooks; we follow Gurtin et al. (2010). 
The principal invariants of the left Cauchy-Green strain tensor B are 

/i=trB, /2(B) = ^[(trB)2-tr(B2)], /3(B) = detB. (36) 

The most general form of the constitutive response function for the Cauchy stress of an incom- 
pressible homogeneous isotropic material in an undistorted reference configuration is 

cr = -(7l + /3i(XB)B + /32(XB)B-\ (37) 

where Ib is the list of only the first two invariants of B (incompressibility demands that I3 = 
1). Let {xi,X2,X3) denote the Cartesian coordinates of a point in the deformed configuration 
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whose Cartesian coordinates in tlic reference configuration are {Xi, X2, X3). A simple shear is a 
homogeneous deformation of the form 

Xi=Xi+(:X2, X2 = X2, X3=X3, (38) 

where ^ = tan ^ > is a prescribed scalar parameter measuring the tangent of the shear angle 6. 
Consider a body that, in its reference configuration, is a rectangular prism with sides a, b, and c 
lying along the coordinate axes Xi, X2, and X3 respectively (one vertex of the prism coincides with 
the origin). For the particular deformation in Eq. (38), in the deformed configuration, (vr/2) — 9 is 
the angle between sides a and 6 and the XiX2-plane is the shear plane. Equations (38) imply that 
the principal invariants of B for the deformation at hand are 

/i(B) = 3 + C', /2(B) = 3 + C', and 13 = 1, (39) 

where the last of Eqs. (39) implies that a simple shear is an isochoric deformation and therefore 
consistent with incompressibility. We note that the stretch in the X3 direction is equal to 1. The 
components of the Cauchy stress of an incompressible homogeneous isotropic material subject to 
the deformation in Eqs. (38), are 
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where /?i and /32 are the expressions for /3i and (32 corresponding to invariants in Eqs. (39). We 
observe that the shear components of stress are unaffected by the Lagrange multiplier q and that 
the '12' components of the tensors B and B~^ are ( and — C, respectively, that is they are a direct 
measure of the shear deformation. The shear modulus is the ratio between the stress component 
associated with the shear plane and the corresponding component of B. Referring to Eq. (38), 
recalling that the XiX2-plane is the shear plane, and that B12 = C, the shear modulus is given by 

/.= f' =/3i(0-^2(0. (41) 

To distinguish this shear modulus determined via a simple shear test from the effective shear 
modulus introduced later, we will refer to n as the "true" shear modulus. 



5.2 Shear response via the inflation and extension test 

For isotropic materials, and in fact even for orthotropic materials if the directions of orthotropy are 
appropriately aligned, the principal directions of strain coincide with the corresponding principal 
directions of the Cauchy stress. In the extension and inflation test, the principal directions of 
deformation are parallel to the coordinate lines of the coordinate system (Fig. 1). Using elementary 
tensor transformations, we deduce that, at each point, the maximum shear strain occurs along lines 
in the principal planes that bisect the coordinate directions. 

Consider a cross section perpendicular to the axis of the cylinder, as shown in Fig. 9, where 
the (^-dependent) lines ^ and rj are defined such that ^ forms a 45° angle with the r axis. The 
orthonormal unit vectors ii^ and tt^ in the ^ and rj directions can be expressed in terms of the unit 
vectors Ur and uq as follows: 

= -^{ur + iie) and = -^{—Uj. + iie). (42) 
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Figure 9: Directions of maximum shear. 



Since B and a have the forms B 



Brr Ur (^ilr + Bog ilg ilg + Bzz Uz ® Uz and CT = Grr ilr ® Ur + 



Ue + (Tzz 



(8) iiz, the ^r] (shear) components of these tensors are 
ii^ ■ Bii^ = ^{Bgg - Brr) and a^n = u^- (^Uri = \{aee 



Or 



(43) 



Note that the shear stress a^^ is completely unaffected by incompressibility and it can be argued 
that it is also unaffected by the presence of a residual stress state.* 

We now introduce the logarithmic average through the wall of the vessel, which we denote by 
\»\ and define as follows: 

1 r° dr 

Using the above definition, along with the last of Eqs. (43) and the equilibrium equation in Eq. (12), 
which is valid for any constitutive model for which the principal directions of stress are parallel to 
the principal directions of strain in the infiation and extension test, we have 



(45) 



21n(ro/ri)' 

where we have used the pressure boundary conditions on the inner and outer surfaces of the cylinder 
(with the reference state taken as the no load state). 

The stress estimate in Eq. (45) is as easy to obtain experimentally as that of current methods 
based on Eq. (3) or its correct form in Eq. (2). However, the estimate in Eq. (45) is completely 
unaffected by incompressibility. In addition, and most importantly, using Eq. (45) we can obtain an 
effective (i.e., not pointwise) yet rigorous measure of a shear modulus of the material, again immune 
from incompressibility effects. Mimicking the definition of shear modulus from the analysis of the 
simple shear test, and denoting the vessel's effective shear modulus by /Ueg, we define ^eflf as follows: 



(46) 



That is, we define effective the (wall averaged) shear modulus as the ratio of the Cauchy shear stress 
to the corresponding Cauchy-Green shear strain, where both are log-averaged over the esophageal 
wall with the lumen represented as an effective cylinder. Using Eq. (10) and the first of Eqs. (43) 
we have that, in an inflation and extension test, B^^ = (l/2)[(r^/i?^) — i?^/(A^r^)], so that 



1 



21n(ro/ri 



r 



W \ dr 



A2r2 



(47) 



*This is certainly the case for isotropic materials in an undistorted configuration. It is also the case for orthotropic 
materials in an undistorted configuration whose directions of orthotropy are parallel to the principal directions of 
strain in the inflation and extension test. 
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As it turns out, the integral in Eq. (47) coincides with the function 7(ri, ro, Ri, Ro, A^) in Eq. (17). 
Therefore, 

so that, substituting Eqs. (45) and (48) into Eq. (46), and using the definition in Eq. (17), 
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1 rRi Rl\^]^JnRo 



(49) 



where, for the class of deformations being considered here and from Eq. (7), A^ = {R^— Rf) / (rg—rf) 
is the ratio of cross sectional wall areas in the reference state relative to the deformed state. For 
isotropic incompressible tissue, this ratio quantifies the stretch in longitudinal direction of the 
esophageal wall relative to the (no-load) reference state,* and A^ < 1 implies longitudinal shortening. 

The effective shear modulus HeS can be used in experiments as a rigorous measure of the elastic 
shear response of any incompressible isotropic material. Since the formula for ^gff is explicitly 
given, its computation does not require the determination of a slope from a graph of stress vs. 
strain. As such, Eq. (49) offers a direct estimation of the elastic response of the material that 
is more straightforward to obtain than that of current approaches. Clearly, |<t^,j1 can always be 
plotted against IB^,^} if desired. However, because the stiffness measure is the ratio of stress to 
strain, it must be interpreted as a secant modulus rather than as a tangent modulus, which is the 
slope of a stress vs. strain curve. 

Remark 1 (Comparison with the simple shear test). The quantity fieS is defined as the ratio of two 
average quantities. Therefore, it pertains to a structure as a whole rather than being the outcome 
of a homogeneous deformation (like the simple shear test), which can be said to apply pointwise. 
With this in mind, to compare fj,es with the true shear modulus, referring to the discussion following 
Eq. (40), we recall that, in the simple shear test, the components B12 and such that 

B12 = —Bi2- Furthermore, the principal stretch in the direction perpendicular to the shear plane 
is equal to 1. A similar situation occurs in the inflation and extension test when A^ = 1, that is, 
{B^^l = —{B^^}. Therefore, we will compare the true shear modulus of a variety of models with 
the correspondent fx^s only for A^ = 1. However, the definition of /leS is not constrained by A^ = 1 
in that it always represents a rigorous measure of the elastic shear response of the material. 

5.3 Comparison between //eff and the true sheeir modulus for various models 

Using the definition of and Eq. (37), we derive expressions for fi^s for some common constitutive 
models. To facilitate these derivations, wc first derive expressions for B and B~^ in terms of the 
principal stretches A^ = R/{Xzr) and (A^ = (A^A,.)^^). From Eq. (10), 

E = y^ihr^ilr + ilQ ® Uq + A^ Uz ® Uz (50) 

A, A„ 



and 



B =—rUr®Ur + ^zK'^e®Ue+-^Uz®Uz- (51) 
A„ A^ 



The corresponding expressions for B^^^ and B^^ are 



*The same could be said for orthotropic materials such that the longitudinal axis of the esophagus is also a 

direction of material symmetry 
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which imply 



B 



Using the definition of /Xefr and Eq. (37) , 



(53) 



(54) 



Neo-Hookean and Mooney-Rivlin Models. Since f3i = ji and /32 = for the neo-Hookean 
model, and = /ii and ^2 = — A*2 for the (full) Mooney-Rivlin model, Eq. (54) yields 



(/^eff)NH = and (/Xeff)MR = /^l + /U2 



(A2^-A2J + (2/A,)ln(A,jA,J- 



(55) 



Since ^ is the true shear modulus of a neo-Hookean material, the first of Eqs. (55) implies that 
(/^eff)NH is exactly equal to the true shear modulus. To compare (/Ueff)MR to the true shear modulus 
of the MR model we let A^ = 1 in the second of Eqs. (55). This gives 



(/^cff)MR| 



(56) 



which, again, is the true shear modulus for the (full) Mooney-Rivlin model. 

Gent Model. We now consider the Gent model (see, e.g., Gurtin et al., 2010): 

-1 



-q\ + 2ci 1 



/1 



Tin 
^1 



B-2c2lo^B-\ 



(57) 



where Ii and I2 are defined in Eqs. (36), and where ci > 0, C2 > 0, and I™ > 3 are constants. The 
true shear modulus of the Gent material is (Gurtin et al., 2010) 



= 2ci( 



Tm 
-'1 



3+e' 

The functions /3i and /32 for the Gent model in an inflation and extension test are 

2ci/i™A2a2 ^ 2C2\\\1 



(/rn + 3 _ ^2)^2^2 _ ^2^4 _ l 

Therefore, the /Lteff for the Gent model is 



and 1^2 = — 



(58) 



(59) 
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l/^ 2AtA^.+l 



^^IA|T+A^7+Aj 



[A2^-A2 +(2/A,)ln(A,„/A,J] 



(60) 



where k = /™ + 3 — A^. To compare (Mefl[)Gent with the true shear modulus in Eq. (58), let Aj; = 1 
in Eq. (60) to obtain 



ilJ-es) 



Gent 



Xz = l 
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^4 - (2 -h /r)2 
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V^4 - (2 + /f )2 
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^4- (2+/f)2y. 
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'^l-(2 + /r)A2 +A^^ 

1 /2A2 +1 
V V3 

+ C2 In \ [Xl - \l + 2 ln(A,„/A,J] (61) 

We compare Eqs. (58) and (61) in Fig. 10, where we have set C, in Eq. (58) equal to [-B^rjl];^ =i in 
the first of Eqs. (49), i.e., 

^ = AHvo/n) - -^'^ + ^ HKJK,)] . (62) 

We non-dimensionalized the expressions for the shear moduli with ci , and have chosen the following 
values of the constitutive parameters: 

— = ^1,2 and /r = 500. (63) 

Cl 

In the Gent model, the quantity /{" serves as an upper bound of the quantity Ii. If /{" is too small, 
it is not possible to explore the behavior of the material in the thin wall limit as we have done 
earlier. Recalling that the physiological range in the esophagus had 0.12 < t/vi < 0.69, we can see 
that the effective shear modulus tracks the true shear modulus well in the physiological range. The 
difference between the two shear moduli is well below 10% even around t/rj = 0.1. To put things 
in perspective, the value t/ri = 0.1 corresponds to iB^^jJ = 2.664, which, in turn, corresponds to a 
shear angle of 69.43°. This level of shear is extreme for the biological systems of interest. 

Fung Model. We now consider an isotropic incompressible Fung-like model of the type (Fung, 
1993) 

o- = -Ql + 6e"(-^i-3)B, (64) 
where a and b are constants and Ii = tr B. The true shear modulus for this model is 

H = he"-^\ (65) 

In an inflation and extension test, the functions /3i and /32 for this model are 

^^ = he<^r+K''K^+xl-^) and /^a = 0. (66) 

As simple as the expression for /3i may seem, we were unable to obtain a closed-form expression 
for the /v,cfr of the Fung model. Therefore, our results are based on a numerical evaluation of 
carried out using Mathematica 8 (Wolfram Research, Inc., 2010). In Fig. 11, we plot n/b from 
Eq. (65) along with i^eS obtained by substituting Eqs. (66) into Eq. (46) for the case with A^ = 1. 
We considered two cases with a = 0.1 and 1. As was done in Fig. 10 for the Gent model, for the 
purpose of comparison, we let ( = [[i?5r;l||;^ -i in Eq. (65). Figure 11 shows that the approximation 
of /X by /Xeff is very good over the entire physiological range. Not surprisingly, the quality of the 
approximation decreases for very small values of t/ri corresponding to extremely large the shear 
strains that are generally outside the physiological ranges of interest. 
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Figure 10: Comparison between the nondimensional true shear modulus and the effective shear 
modulus for the Gent Model. The nondimensionalization was carried out relative to the c\ param- 
eter appearing in the Gent model. The solid blue lines correspond to Eq. (58) with C, = \B^ril and 
are denoted by ji. The solid red lines correspond to Eq. (61) and are denoted by //eff- The two cyan 
lines delimit the physiological range. 



6 Example: a re-evaluation of the standard method from the lit- 
erature 

We have shown that the slope of total hoop stress plotted against strain does not approximate a 
true material stiffness in principle and that the Laplace law does not approximate hoop stress for 
typical tubular vessels in the human body in practice. To illustrate our points, we here re-evaluate 
data given in Takeda et al. (2003) using our proposed method to estimate a true effective shear 
modulus with the same data collected for the Laplace law method. We used Takeda et al. to 
illustrate our method primarily because they have presented data in plots and statistical measures 
that make possible a re-evaluation and comparison using our proposed method. 

We refer to the Laplace-law-based approach applied by Takeda et al. as the "standard approach." 
The aim of the experiments by Takeda et al. was to determine the effects of atropine, a smooth 
muscle relaxant, on esophageal wall stiffness, where "stiffness" was quantified by applying the 
standard approach to estimate what was thought to be an extensional modulus of elasticity (see 
Section 1), as described here below. 

In Takeda et al. (2003), intraluminal pressure was measured in the esophagus with a mano- 
metric catheter at the same location as an intraluminal ultrasound probe that imaged the cross 
section of the esophageal wall, ~ 7.5 cm above the lower esophageal sphincter (high pressure zone). 
Approximating the pressure external to the esophageal wall as atmospheric and estimating the 
thickness of the esophageal wall from the ultrasound images by mapping the cross sectional areas 
into circles with equivalent areas, the average total circumferential stress of the esophageal wall 
was estimated using the Laplace equation, namely Eq. (3). To estimate wall stiffness, Takeda et al. 
collected data with systematic distension of the esophageal lumen by progressively filling a high- 
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Figure 11: Comparison between the nondimensional shear modulus and the effective shear mod- 
ulus for the Fung model. The nondimensionalization was carried out relative to the b parameter 
appearing in the Fung model. The solid blue lines correspond to Eq. (65) with = [[-B^?y]] and are 
denoted by fi. The solid red lines correspond to Eq. (54) for the Fung model and are denoted by 
Hes- The two cyan lines delimit the physiological range. 



compliance bag within the lumen with water. The bag was attached to the manometric catheter so 
that bag pressure and volume was measured simultaneously during distension. Using the changes 
in intraluminal pressure from manometry and changes in cross sectional geometry estimated from 
ultrasound images, wall hoop stress was quantified using the Laplace law and circumferential strain 
was quantified, relative to the resting state, as relative change in effective circumferential dimension 
of the esophageal lumen. Takeda et al. then used slopes of stress against strain plots as an elastic 
modulus to compare esophageal wall properties before vs. after administering atropine. 

Takeda et al. collected their distension data using two different protocols. In one protocol, the 
bag was filled in steps by systematically increasing bag volume (isovolumic distension) from 5 to 
20 ml. In the second protocol the bag was filled by systemically increasing bag pressure (isobaric 
distension) from 10 to 60 mmHg. The level of distension and the time over which the lumen was 
distended in each step were different for the two protocols. The rate of distension and the time 
delay before data were collected, potentially affecting the generation of active tone, were not given. 
Nevertheless, one expects similar changes in wall properties in response to atropine, independent 
of protocol. 

For 8-10 asymptomatic subjects the following quantifications were estimated at each distension 
by Takeda et al.: intraluminal pressure, lumen cross sectional area, outer wall circumference, and 
approximate wall thickness. The ratio of outer circumference in the distended vs. resting state 
was computed at each distention and was used to convert the Laplace law hoop stress to the 
corresponding component of the second Piola-Kirchhoff stress and to calculate the corresponding 
component of the Lagrangian strain. Plots of (aforementioned components of) the second Piola- 
Kirchhoff stress against Lagrangian strain were fit with straight lines over the ranges of applied 
distension, the slopes of which were interpreted as moduli of elasticity. 
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A primary conclusion from Takeda et al. was stated in the abstract: 

The Youngs modulus, which is the slope of a linear stress-strain relationship, was sig- 
nificantly higher after atropine in the isovolumic study but not in the isobaric study. 

Specifically, Takeda et al. estimated the Youngs modulus to be 37mmHg vs. lOTmmHg before vs. 
after atropine with isovolumic distension, while with isobaric distension they estimated 102mmHg 
before atropine vs. lOTmmHg after atropine was applied, a change that was not statistically sig- 
nificant. This result lead to the following conclusion (see the discussion section of Takeda et al., 
2003): 

These results suggest that atropine behaves as if it stiffens the esophageal body in an 
isovolumic but not an isobaric study, 

followed the obvious question 

Why does atropine behave in different ways depending on the modality of esophageal 

distension? 

Takeda ct al. go on to argue that their result implies a protocol-dependent active response to 
atropine overlying a protocol-independent passive response. 

However, we have argued in this paper that total stress response to strain is seriously polluted by 
incompressibility in the hydrostatic component, that the tangent to the curve of a total principal 
stress component against its corresponding strain component therefore does not provide a true 
elastic modulus, and that the standard approach is fundamentally erroneous, so that a true modulus 
requires filtering out the hydrostatic component from total stress. In addition, we have shown that 
Laplace law estimates for total stress are seriously in error for the esophagus and other thick- 
walled elastic biological vessels, and that the Laplace law stress cannot be used as a surrogate 
for the strain-dependent component of total stress. We re-evaluate the Takeda et al. result with 
our approach based on an effective shear stress that is unpolluted by incompressibility and avoids 
the Laplace law. With this re-evaluation, we also point out differences in application between our 
proposed method and the standard approach. 

6.1 Application of the new method to the data by Takeda et al. (2003) 

Whereas the standard approach uses the slope of total hoop stress against strain to define an 
elasticity modulus, our approach estimates an effective shear modulus directly with Eq. (49). Like 
the standard approach, we replace the true lumen cross section by an effective circular cross section 
with equivalent cross sectional areas. The outer and inner radii of the effective lumen are determined 
in the reference configuration (i?j,i?o) and under distension {ri,ro)- To apply Eq. (49) one first 
quantifies the reference configuration, assumed to be the no-load state (i.e., the state with zero 
transmural pressure difference, Ap). Then, with an inflation experiment (isobaric or isovolumic), 
the transmural pressure difference Ap is quantified together with inner and outer radius rj and 
of the effective lumen at each distension. Since manometry measures only intraluminal pressure, 
to estimate Ap the external pressure is assumed to be close to atmospheric. Finally, longitudinal 
stretch in the wall material must be estimated from the second of Eqs. (7), which is a statement 
of incompressibility during distension from {Ri,Ro) to (ri,ro). With these quantifications, one 
applies Eq. (49) to estimate effective shear modulus. 

Note that the data required to estimate and effective shear modulus are the same as the data 
required for the standard approach: inner and outer radii needed in Eq. (49) are also needed to 
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quantify hoop stretch in the form of rj and wall thickness (r^ — rj). With our new method, however, 
particular care should be taken to quantify the geometry of the reference state as accurately as 
possible. 

Whereas the same information needed in the standard approach is also applied in our new 
method, our new method requires the user to quantify explicitly the local longitudinal stretch Aj, 
a physiologically important quantity worth studying in its own right (see, e.g., Mittal et al., 2006; 
Nicosia et al., 2001; Shi et al., 2002), and to incorporate the quantification of longitudinal stretch 
into our effective shear modulus. The fact that the traditional approach does not use longitudinal 
stretch to estimate stiffness is a manifestation of not removing the important confounding effects of 
incompressibility from total stress to produce the stress component that represents clastic response. 

To extract the information we needed for the application of our method, i.e., to apply Eq. (49), 
we digitized data presented in plot form in Takeda et al. (2003). It should be recognized, however, 
that had the experiment been designed to directly apply Eq. (49), it would have been developed 
and analyzed somewhat differently, and we would have had available the raw pre-analyzcd data. To 
obtain the necessary quantifications, we digitized the ensemble- averaged results presented in Figs. 1, 
2, and 3 of Takeda et al. In the isovolumic case, we quantified pressure differences from Fig. lA, 
assuming that Takeda et al. approximated intrathoracic pressure external to the esophageal wall 
as atmospheric, as is typical. From Fig. IB we quantified average intraluminal cross sectional area 
(CSA). 

Although the esophageal lumen cross section is not circular in general (cf. Nicosia et al., 2001), 
to estimate effective properties of the lumen wall, we placed the CSA within a circular cross section 
to quantify inner radius so that rf = vr/CSA. We then used the data from Fig. 3A to obtain 
measurements of the hoop component of the second Piola-Kirchhoff stress a (according to the 
notation in Takeda et al., 2003) and the corresponding Lagrangian strain e (Eq. (3) in Takeda et al., 
2003). From e, we obtained the corresponding circumferential stretch ratio, using e = |(Ag — 1), 
as given in Eq. (3) of Takeda et al. From the Laplace law relationship between the hoop component 
of the second Piola-Kirchoff stress and wall thickness used by Takeda et al., specifically their Eq. (5), 
we were able to extract the wall thickness h and, from that, rg = ri + h. To estimate the outer 
radius in the reference state Ro, we used the averages for outer circumference supplied by Takeda 
et al. (see p. 624) before and after the administration of atropine. To obtain Ri, we applied a linear 
regression to the thickness vs. pressure data and extrapolated to Ap = 0. Combining these results 
with Ro yielded the desired inner radii values in the reference configuration. A similar strategy was 
used to analyze the data in the isobaric case. 

Since the analysis developed in this paper was based on the averages presented in Takeda et al. 
(2003), we carried out a sensitivity analysis within the range of values provided by Takeda et al. 
for the outer circumference in the rest state. We found that the essential results are insensitive to 
the uncertainly in circumference and the conclusions are unchanged. 

6.2 Analysis of the Effective modulus with the new method 

As discussed above, to quantify the effective shear modulus via Eq. (49), it is necessary to estimate 
the local longitudinal shortening of the muscle layers given by Xz in Eq. (7) evaluated at the 
outer surface of the lumen. The quantity A^ measures the local longitudinal dimension of the 
lumen relative to the longitudinal dimension in the reference (no-load) state. Relative shortening 
(lengthening) of that axial location in the lumen is quantified by reduction (increase) in A^ from 
1. Using the data from the plots in Takeda et al. (2003) as described above, the variation of Xz 
due to isobaric (pressure control) and isovolumic (volume control) distension of the water-filled 
bag, as described above, is shown in Fig. 12, in which Xz is plotted as a function of intraluminal 



23 



pressure relative to atmospheric pressure. Figure 12 indicates that the esophageal lumen at the 
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Figure 12: Local longitudinal shortening (A^) of the esophageal wall plotted against transmural 
pressure difference Ap for isovolumic and isobaric lumen distension before (BA) and after (AA) 
atropine was administered, as per protocols and data presented in Takeda et al. (2003). 

measurement location did shorten longitudinally in response to distension, indicating the activation 
of longitudinal muscle fibers. The degree of shortening increased with increasing distension at the 
same rate with both protocols, both before and after atropine was administered. The isobaric 
protocol, in which pressure in the bag was increased stepwise, suggests that the level of longitudinal 
muscle shortening is about the same before and after atropine, while the isovolumic distension result 
is less clear since the relationship between pressure and distension was significantly altered by the 
smooth-muscle-relaxing effects of atropine (Takeda et al., 2003). This is indicated in Fig. 12 by 
the shift in the post-atropine curve to the left relative to the pre-atropine curve. The results 
also seem to suggest that the level of longitudinal shortening at the same intraluminal pressure 
is protocol dependent, with greater longitudinal shortening under pressure control than volume 
control. However, the relative shifts in the curves change when Xz is plotted against distention (r,) 
rather than pressure, so one cannot draw general conclusions from the comparison of magnitude of 
longitudinal shortening between protocols and before vs. after atropine. Nevertheless, distension 
stimulates longitudinal shortening in all cases and the level of shortening increases with increasing 
distention at about the same rate for both protocols and before vs. after administering atropine. 
We conclude, therefore, that the stimulation of longitudinal shortening at the same rate before vs. 
after administering atropine is a general result, independent of protocol. 

The effective shear modulus of the esophageal wall, Hes, is plotted in Fig. 13 against the effective 
shear component of the left Cauchy-Green strain, Egg = IB^^J, using the data in Takeda et al. 
(2003) and the results for longitudinal shortening parameter in Fig. 12. Comparisons are made 
before vs. after the administration of atropine for the isobaric and isovolumic protocols in which 
bag pressure or bag volume were controlled during distension, and comparing effective wall stiffness. 
Note that the two after atropine (AA) results become unstable in the limit of low strain, when 
-Beff ~ 0.5. This is because in Eq. (49) both the numerator and the denominator approach zero 
in the reference state limit, since Egg and Ap are both zero in the reference state by definition. 
To maintain finite /Ueff as i?efr 0, the denominator must approach zero at the same rate as the 
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Figure 13: Effective shear modulus of the esophageal wall /Xeg (Eq. (49)) plotted against the com- 
ponent -Beff = I-B^r^l of the effective left Cauchy-Green strain in Eq. (48) for (a) isovolumic and (b) 
isobaric lumen distension before (BA) and after (AA) atropine was administered. 



numerator; however the ratio is sensitive to imprecision in the measured data and truncation error 
in the no-load limit. The unstable range must therefore be excluded in the analysis of ^eff; and in 
Fig. 13 we only analyze the range Seff > 0.5. 

The standard method used by Takeda et al., where an apparent modulus is deduced by plot- 
ting the Laplace law approximation of total hoop stress against strain, led to inconsistent results 
depending on which protocol was used to distend the esophageal wall: the apparent modulus in- 
creased after atropine when the esophagus was distended by incremental increases in bag volume, 
but did not change when the esophagus was distended by increasing bag pressure. The apparent 
result that the effect of atropine on esophageal wall stiffness is protocol-dependent led Takeda et al. 
to suggest that the tonic response of esophageal muscle is different depending on how the esophagus 
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is distended. Atropine is a smooth muscle relaxant; it is unclear why the change in muscle stiffness 
from unrelaxed to relaxed state would depend on the method used to distend muscle. 

Figure 13, however, shows that the inconsistency disappears when the standard approach is 
replaced by the application of the effective shear modulus developed in Section 5.2: Fig. 13 shows 
that the effective shear modulus increased after administering atropine in both the isovolumnic and 
isobaric protocols. The change in muscle stiffness with atropine is independent of the protocol used 
to distend the esophageal wall. We conclude that the resistance to distension is greater when the 
esophageal smooth muscle is in an atropine-induced relaxed state than in the normal state in which 
the response is a summation of passive elastic resistance and changes in active tone stimulated by 
stretch. This conclusion, based on a true measure of elastic modulus, is a general one, independent 
of protocol, resolving the previously inconsistent result. The resolution of inconsistency, we argue, 
is a consequence of removing the two essential errors in the standard method: (1) the pollution of 
the hydrostatic component of the total stress response to strain by incompressibility, and (2) the 
inaccurate representation of both total and clastic stress by the Laplace law. 

Having generalized the essential result that esophageal muscle is stiffer after atropine, we can 
note potential application-dependent differences in the results of Fig. 13. In particular, although 
Takeda et al. using the standard method concluded that the elastic modulus was independent 
of the magnitude of distension in all cases. Fig. 13 suggests that this may be the case for the 
isovolumic protocol, but that with the isobaric protocol the esophageal wall appeared to stiffen 
with increasing distension, both before and after atropine was administered. Additionally, Takeda 
et al. report apparent elastic moduli 37-83 mmHg before-after atropine with isovolumic distension 
and ~ 105 mmHg with isobaric distension, while our shear modulus estimates are ~ 50-90 mmHg 
before atropine and « 100-150 mmHg after atropine. Whereas the standard method indicated 
an increase in elastic modulus of 46 mmHg after atropine with the isobaric protocol and no 
change with isovolumic distension, we estimate a change in effective shear modulus of 30-50 mmHg 
depending on distension level and, perhaps slightly, distension protocol. However without access 
to the raw data and a more careful statistical analysis, one should not draw firm conclusions from 
these details. 

Our overall conclusions are (1) that the standard method does not provide a proper measure 
for elastic modulus due to pollution of the hydrostatic component of the total stress due to in- 
compressibility and serious inaccuracies in the Laplace law approximation of total hoop stress, and 
(2) that our proposed method based on shear stress alleviates the problem by isolating the elastic 
component of stress and providing a method for estimating stress without use of the Laplace law. 

7 Summary and Conclusions 

In this paper we discuss the applicability of a current common approach to estimate total average 
hoop stress across the wall of biological vessels, and the practice of estimating material stiffness 
as a modulus of elasticity, from a plot of total hoop stress against circumferential strain. We 
refer to this combination as "the standard method." We have shown that both elements of the 
standard method, the use of total stress estimate an elastic modulus when plotted against strain 
and the application of the Laplace law approximation to estimate total hoop stress, are inapplicable 
for biological vessels which are relatively thick and are composed of incompressible tissue such as 
muscle. 

The Laplace law provides an incorrect estimate of the transmural average of the total hoop 
stress because the crucial assumption, that the relative thickness of the vessel wall is very small 
in comparison to the relative transmural pressure difference, is rarely met in biological vessels. 



26 



In comparison to bubbles, for which the Laplace law was originally intended, the Laplace law 
does not apply to biological vessels as a consequence of the relatively thick nature of the vessel 
walls, and the deformability of biological tissue in response to relatively low transmural pressure 
differences. In contrast, the Laplace law is generally applicable to finite-thickness metal tubes when 
they arc supporting very large relative pressure differences. However, in any application for which 
an estimate of the transmural average of total hoop stress is sought in a human or animal vessels, 
we argue in this paper that the Laplace law should be avoided. 

The Laplace law is a straightforward reduction of a simple equilibrium force balance where 
the pressure external to the vessel appears only within the transmural pressure difference and 
the independent effect of external pressure is neglected. However, it is common to observe the 
Laplace law as a force balance reduction with the external pressure set to zero in all terms. This 
approximation is not valid in biological vessels and the common practice of setting the external 
pressure to zero in the Laplace law should be avoided. The Laplace law estimate of the average 
total hoop stress is easily corrected by simply including the external pressure in the force balance 
and not making the "thin wall approximation." The common practice of estimating the elastic 
modulus of the wall from the slope of the average total hoop stress plotted against an appropriate 
strain measure should remain valid since, from Eq. (2), the thin wall approximation simply shifts the 
average stress response by the amount po, thus leaving unaffected the slope of this response function. 
However there is a more serious issue in this practice that invalidates the interpretation of slope as 
a response to the elastic deformation: the incompressibility of muscle contaminates the total stress 
with a contribution that is largely independent of the elastic response of the material. To interpret 
the slope of total stress against strain as elastic modulus, one assumes that the contribution due to 
incompressibility is either independent of deformation or sufficiently small relative to total stress. 
We have shown that incompressibility can be a major contributor to total hoop stress (Fig. 6). 
Thus, in general, a plot of transmural average of the total hoop stress against strain cannot be used 
with confidence to estimate a modulus of elasticity, even in an approximate fashion. 

We have proposed an alternative approach to overcome the difficulties inherent to the use of the 
average total hoop stress to estimate a modulus of elasticity. The proposed approach uses the same 
information typically gathered during an inflation and extension test to derive a mechanistically 
rigorous measure of the shear response of the material. It is a practical method that, like the current 
approach, can be applied in vivo as well as in vitro. The rationale for this approach is that the 
shear response for isotropic materials and, in some cases, for orthotropic materials, is completely 
immune from the effects of incompressibility. 

The proposed approach provides a measure of an effective shear modulus of the material. We 
have shown via the analysis of several common models for incompressible isotropic materials that 
the proposed approach delivers very good estimates of the shear modulus of the material especially 
over the range of deformations that are of interest in the characterization of biological tissues. 

We further re-analyzed data from Takeda et al. (2003) presented in their plots and statistics. 
Takeda et al. used the standard method (deduction of a modulus from the slope of total hoop stress 
plotted against strain and using the Laplace law to approximate total hoop stress) to estimate an 
apparent modulus of elasticity with two different protocols to distend the esophageal lumen, one 
in which the volume of a water filled back in increased in controlled incremental fashion, and one 
in which water is added to the bag to increase the pressure in the bag in controlled increments. 
They found that the change in stiffness after administering atropine, a smooth-muscle relaxant, was 
different depending on the protocol used to distend the esophagus: the esophageal wall apparently 
stiffened after atropine in one protocol but not the other. We showed that this inconsistency is 
resolved when our new method for estimating a true shear modulus is applied: the elastic shear 
modulus is observed to increase after atropine is administered under esophageal distension with 
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either protocol. 

We argue that the previous "standard" method to estimate a modulus of elasticity of biological 
vessel walls based on total hoop stress and the Laplace law should be abandoned. We further suggest 
that our method is a useful replacement for the standard method. Our proposed new approach 
provides a true measure of elastic modulus without the polluting effects of incompressibility in 
methods based on total hoop stress and without application of the Laplace law which is highly 
inaccurate for biological vessels. A further major advantage of our approach is that the effective 
shear modulus of the material can be measured directly from experimentally obtained quantities. 
This is different from current approaches, whereby a modulus is inferred from the slope of an 
experimentally determined stress-strain curve. In this sense, the proposed approach offers a direct 
accurate measure of the elastic shear response of incompressible isotropic materials. 



Appendix 

The formulation of the problem concerning the inflation and extension of an isotropic nonlinear 
elastic circular cylinder can be found in several textbooks (see, e.g., Humphrey, 2010; Ogden, 1997). 
Aspects of the solution to this problem are also available. However, the authors could not find in the 
literature the specific expressions they needed in this paper for the Lagrange multiplier, the Cauchy, 
and second Piola-Kirchhoff stress tensor and their averages through the thickness. Therefore, we 
present our derivation of the quantities in question. We do not claim to have found a new solution 
to the inflation and extension problem. We have simply derived the expressions for the quantities 
needed by our analysis. What is new is our paper is the definition of an effective shear modulus of 
a tubular vessel that can be measured using current experimental approaches both in vivo and in 
vitro and, most importantly, that properly isolates the elastic response of tubular biological vessels 
from the response to incompressibility. 



Inflation and extension of a neo-Hookean tubular vessel 

Recalling that I = u,- (S) u,- + uq $5 uq + Uz'Si iiz and that B is given in Eq. (10), for a new-Hookean 
material described via Eq. (14), we have 

a^^ = -g + — — and aee = -q + IJ-j^- (A.l) 

Substituting Eq. (A.l) into Eq. (12) (the only nontrivial component of the balance of linear mo- 
mentum for the problem at hand), we have 

The deformation under consideration allows us to express r as a function of R and vice versa. When 
we view R = R{r), from the first Eqs. (7), we conclude that dR/dr = X^r/R. This allows us to 
write the last two terms of the left-hand side of Eq. (A. 2) as follows: 

R^ _ _l d R^ RdR _ I d R^ ^ dlnr r _ 1 dlni? 

j,3 2 dr dr 2 dr r^ ^ dr R'^ dr 

Hence, Eq. (A. 2) can be given the form 

dq /i d i?^ /i d , r ^ R^ ^ ^ r /* ,x 

dr = 2^Tr7^^XzTr^''R ^ ^ = ^ + 2AI7^ + i?' ^^'^^ 
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where K is a constant of integration that wc find by enforcing a boundary condition. The outer 
boundary of the cylinder, with outward unit normal iir, is subject to the pressure po- Hence, 
aUr\r=ro = ~PoUr, which implies that arr{ro) = —Po- This expression, along with the first of 
Eqs. (A.l) and the last of Eqs. (A. 4), yields an equation for K whose solution is 



(A.5) 



Substituting Eq. (A.5) into the last of Eqs. (A.4) and simplifying, we have 



+ 



(A.6) 



We now observe that the the inner boundary of the cylinder is subject to the pressure Pi, which 

implies that arr{ri) = —pi- We can use this expression along with Eq. (A.6) to determine an 
expression relating the transmural pressure difference Ap to the elastic constant Doing so gives 
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Recalling that Ap = Pi— Po, the above equation can be solved for /x to obtain 



/X = Ap 



J_( Rl Rj 
2\l \ rl rf 



^ 1 ^JnRo 

Xz \roRi 



(A.8) 



The full Cauchy stress solution is then obtained by substituting the expression for q into Eq. (14) 
and using Eq. (10). This gives 



cr = — 



H {rRo 
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1 E? 
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Using the first of Eqs. (35), recalling that in our problem det F = 1, and using Eqs. (8) and (9), we 
then have that the full solution for the second Piola-Kirchhoff stress tensor is 



S = - 



Po + 
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(A.10) 



Inflation and extension of a Mooney-Rivlin tubular vessel 

Recalling that B"-*^ is given in Eq. (11), for a Mooney-Rivlin material described via Eq. (22) with 
jii = 0, we have 

(j„ = -q-n2\l^ and aee = -q - ■ (A.ll) 

Substituting Eq. (A. 11) into Eq. (12) (the only nontrivial component of the balance of linear 
momentum for the problem at hand), we have 



dq 



d 



0. 



(A.12) 
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Using Eqs. (A. 3), we can rewrite Eq. (A. 12) as follows: 



1 d i?^ 



1 d 
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dr V A2 r2 



q = C + iJ,2X 



I 1 , r 
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(A.13) 



where C is a constant of integration that we find by enforcing a boundary condition. Recalling that 

on the outer boundary of the cylinder we must havearriro) = —Po, from the first of Eqs. (A. 11) 
and the last of Eqs. (A.13), we obtain an equation for C whose solution is 



C = Po + 1^2X2 



2X1 rl^ Xz"" r„ 



(A.14) 



Substituting Eq. (A.14) into Eq. (A.13) and simplifying, we have 



q = Po + A*2A 
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We now observe that the the inner boundary of the cylinder is subject to the pressure pi, which 
implies that arr{ri) = —pi- We can use this expression along with Eq. (A.15) to determine an 
expression relating the transmural pressure difference Ap to the elastic constant fi2 ■ Doing so gives 
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Recalling that Ap = Pi — Po, the above equation can be solved for 112 to obtain 
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The full Cauchy stress solution is then obtained by substituting the expression for q into Eq. (22) 
and using Eq. (11). This gives 
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Using the first of Eqs. (35), recalling that in our problem det F = 1, and using Eqs. (8) and (9), we 
then have that the full solution for the second Piola-Kirchhoff stress tensor is 
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Formulcis for the Determination of Averages 



The analysis in the body of the paper uses that average through the vessel's wall of a number 
of quantities. As it turns out, the determination of these averages consists of combinations of 
averages of elementary functions. For the sake of a compact presentation, we limit ourselves to the 
elementary functions in question and their corresponding averages. 

Let p = \Jrf — B^/Xz- Recall that B? = X^r^ — X^rf + i?f , which can be rewritten as = 

Xz{r'^ — p^)- The first two elementary functions we consider are 



= AJ 1- 



Their averages are given by, respectively. 
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Concerning logctrithmic averages 

There are a number of expressions that can be written in terms of the principal stretches Aj., Ag, 
and Xz- In addition, due to the incompressibility constraint, these expressions can be rewritten in 
terms of Xz and either Ar or Xq. When computing logarithmic averages of these expressions, it is 
often convenient to express integrands in terms of the principal stretch Xq (and A^, although A^; 
appears as a constant in these averages). In this case, a logarithmic average can be computed as 
follows: 

1 r° dr 1 f'^"" 1 

H = , , , , / • - = , , , , / • vn rT2TdAe, (A.26) 

hi{ro/ri) J^. r ln{ro/ri) J^g. Xo{l - XzXj) 

where 

X0, = ^ and Xe„ = ^. (A.27) 
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